Effect of small particles on the near- wall dynamics of a large particle in a highly 

bidisperse colloidal solution 
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We consider the hydrodynamic effect of small particles on the dynamics of a much larger particle 
moving normal to a planar wall in a highly bidisperse dilute colloidal suspension of spheres. The 
gap ha between the large particle and the wall is assumed to be comparable to the diameter 2a of 
the smaller particles so there is a length-scale separation between the gap width ho and the radius 
of the large particle b ;§> ho- We use this length-scale separation to develop a new lubrication theory 
which takes into account the presence of the smaller particles in the space between the larger particle 
and the wall. The hydrodynamic effect of the small particles on the motion of the large particle is 
characterized by the short time (or high frequency) resistance coefficient. We find that for small 
particle-wall separations ho, the resistance coefficient tends to the asymptotic value corresponding 
to the large particle moving in a clear suspending fluid. For ho ^> a, the resistance coefficient 
approaches the lubrication value corresponding to a particle moving in a fluid with the effective 
viscosity given by the Einstein formula. 



o 

■4-* 



i 

C 

o 
o 



> 

00 

o 

O 

(N 
O 
OO 
O 



X 



I. INTRODUCTION 

Polydisperse suspensions have a significant presence in 
many important industrial and natural systems. Exam- 
ples include dairy products, polymeric gels and nearly all 
complex fluids in biological organisms. Since the dynam- 
ics of colloidal systems is often influenced by their con- 
fining boundaries [l], Q , the confining effects are impor- 
tant in many emerging applications, e.gj in self-assembly 
technology involving colloidal crystals [J, 0] and in chem- 
ical deposition processes on patterned wall Hence, it 
is necessary to understand how polydisperse suspensions 
behave in the proximity of a bounding surface. 

The equilibrium state of confined systems can be fully 
determined by the particle-wall and interparticle interac- 
tion potentials, whereas a description of non-equilibrium 
phenomena is much more complex. In this paper we 
study the non-equilibrium dynamics of a wall-bounded 
bidisperse solution where a large colloidal particle moves 
in the direction normal to the wall in the presence of 
smaller particles. Our objective is to find the correc- 
tion to the near-contact hydrodynamic force on the large 
particle due to the effect of the smaller particles in the 
low-concentration regime. 

One way to tackle this problem is to treat the suspen- 
sion in the gap between the wall and the large particle as 
a fluid with the effective viscosity given by the Einstein's 
formula 6]. This approach, however, is valid only if the 
dimension of the small particles is much smaller than the 
characteristic flow length scale which is the minimum dis- 
tance of separation between the surfaces of the wall and 
the large particle. When the separating distance is com- 
parable to the dimension of the smaller species, a more 
detailed analysis is necessary to obtain the proper cor- 
rection to the hydrodynamic friction. Our main purpose 
is to evaluate this correction in the short time (or high 
frequency) regime, where the influence of the small par- 



ticles on the large particle motion results entirely from 
hydrodynamic effects. 

This article is organized in the following way. In sec- 
tion HH we define our system and discuss assumptions 
involved in our calculations. In section Hill we develop a 
lubrication theory for near-contact approach of the large 
particle towards the wall and use this theory to derive 
expressions for modified lubrication pressure and hydro- 
dynamic force on the large particle in the presence of 
small particles. In section IIV1 we express the dipolar 
source density in the lubrication equations in terms of 
the hydrodynamic forces induced on the surfaces of small 
particles. In section [V] the results of our theory are used 
to determine the resistance force acting on the large par- 
ticle. Our conclusions are drawn in section IVIl 



II. 



HIGHLY ASYMMETRIC BIDISPERSE 
COLLOIDAL SUSPENSION 



We consider a bidisperse colloidal suspension near a 
solid planar surface. Particles of both species are as- 
sumed to be spherical, with high asymmetry in size. We 
describe the influence of the smaller particles on the mo- 
tion of a single larger particle in the direction normal 
to the wall. Our analysis is focused on the effect of the 
smaller species when the separation between the larger 
particle and the wall is comparable to the diameter of 
the smaller particles. 

We concentrate on the short-time (or high-frequency) 
regime; the effect of Brownian motion on the system dy- 
namics can therefore be neglected on the assumption that 
the deviation of the particle distribution from equilibrium 
is small. Under these conditions, the smaller particles af- 
fect the motion of a large particle only by influencing the 
velocity and the pressure fields in the suspending fluid. 
In order to determine these hydrodynamic fields, we first 
analyze the flow in the presence of a single small particle 



FIG. 1: Near-contact configuration of a large particle in the presence of a small particle in a bidisperse colloidal solution. 
Dotted lines represent the boundary of the excluded volumes due to electrostatic repulsion. 



in the gap between the wall and the large particle. The 
cumulative effect of small particles in a dilute suspension 
is obtained by simple statistical averaging over the po- 
sition of the smaller particles. This approach is valid in 
the concentration range where the suspension is dilute 
enough to neglect mutual interactions between smaller 
particles and dense enough for meaningful statistical av- 
eraging. 



coefficient is obtained as the equilibrium average of the 
microscopic stresses produced by the particle motion. 
The results for a general interaction potential can thus 
be obtained by evaluating an appropriate equilibrium av- 
erage of the hydrodynamic functions derived in the fol- 
lowing sections. Our explicit numerical results for the 
excluded-shell potential illustrate the essential features 
of the near-wall dynamics of large particles in a highly 
bidisperse colloidal suspension. 



A. System definition 

The geometry of our system is represented in Fig. [T] 
The small and the large particles have spherical hydro- 
dynamic cores of radii a and b, where o « i. In typi- 
cal highly asymmetric bidisperse colloidal mixtures b/a 
ranges from 10 2 to 10 3 . 

A rigid wall is located at z = 0, and the center of the 
large sphere is on the z-axis. Therefore, the surface of 
this particle can be described by z — h(p) where p is the 
distance from the axis of symmetry (z-axis). The mini- 
mum distance of separation between the particle surface 
and the wall is h(Q) = ha. We consider ho ~ 2a. 

A repulsive screened electrostatic potential between 
the particles is modeled by additive excluded shells sur- 
rounding hydrodynamic core of each particle. The thick- 
nesses of the excluded shells for the small and the large 
spheres are 5 a and <5&, respectively. The thickness of the 
excluded layer associated with the wall is denoted by S w . 
It is assumed that S a ,5b,5 w -C b. 

The interaction potential affects the short-time dy- 
namics of the system only through its influence on the 
equilibrium particle distribution — the effective resistance 



B. Equation for the flow field 

In our calculations, velocity of the large colloidal parti- 
cle is denoted as V p which is in the direction perpendicu- 
lar to the wall. In this paper, all of our results correspond 
to the case where the large particle approaches to the 
wall — the reversal in this motion can be accounted for by 
simple reversal of computed velocity fields and pressure 
gradients. The velocity V p is small so that the creeping 
flow approximation is valid. Hence, the hydrodynamic 
fields are governed by the Stokes equations, 

-VP + ?7V 2 u=0, V-u = 0, (1) 

where u is the velocity field, P is the pressure field, and 
r\ is the viscosity. No-slip boundary condition is assumed 
at the wall and at the surfaces of all particles. 

The hydrodynamic effect of the small particles on the 
surrounding fluid can be described in terms of the in- 
duced force density F on surface of the small spheres. 
Such force density imposes a discontinuity in the gradient 
of the velocity field at the particle-fluid interface. This 
discontinuity is such that we can consider the volume 
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inside the particle as an extension of the fluid domain — 
the fluid in this region undergoes rigid-body motion cor- 
responding to that of the rigid particle. Therefore, the 
entire domain (both inside and outside the particle) can 
be treated as a continuous fluid medium. As a result, 
the flow field described by Eq. |T]) in the presence of 
the smaller particles is equivalent to the flow field which 
obeys the inhomogeneous Stokes equations 

-VP + r/V 2 u = F, V-u = 0. (2) 

Hence, if we can determine F and solve Eq. ^j, we can 
find the effect of the smaller species on the velocity and 
the pressure fields. Accordingly, in the first part of our 
analysis, we solve Eq. to find P and u in the gap 
between the large particle and the wall for an arbitrary 
induced- force distribution F. This analysis is presented 
in section UTT] 

C. Statistically averaged force density 

In order to quantify the cumulative influence of smaller 
species on the hydrodynamic resistance of the large parti- 
cle, we solve an ensemble-averaged flow equation (where 
the average is taken over all possible positions of the small 
particles). This reduced description requires the averag- 
ing of the force density F. Such averaging is valid if 
two conditions are satisfied. Firstly, there should be a 
sufficient number of small particles in the gap between 
the large particle and the wall so that the statistical de- 
scription is physically meaningful. Secondly, the small 
particles should be well separated and mutually non- 
interacting. These two conditions are satisfied for 

ho « C"3 < v/ftfo, (3) 

with C being the number density of the small particles. 
Under these conditions, the statistically averaged force 
distribution can be evaluated from the cumulative effect 
of mutually noninteracting individual small particles. 

III. LUBRICATION SOLUTION IN THE 
NEAR-CONTACT REGION 

In our hydrodynamic calculations, we exploit the fact 
that for near-contact motion of the large particle the 
characteristic length scale along the wall is much greater 
than the characteristic length scale along the perpendic- 
ular direction z. As a result of this separation of length 
scales, we can apply lubrication theory to determine the 
flow and the pressure fields. 

A. Lubrication expansion 

The predominant hydrodynamic effect on the large 
particle near a wall comes from the flow field in the near- 
contact lubrication region. In this lubrication zone, the 



flow can be locally approximated as a flow between two 
rigid parallel planes. The extent of the lubrication re- 
gion is defined by the lubrication length I, which is the 
characteristic length in the direction along the wall: 

/ = V2ab ~ > ho ~ 2a. (4) 

The corresponding characteristic length for the direction 
perpendicular to the wall is the slowly varying gap width 
h ~ 2a <C I. The small parameter that characterizes 
length-scale separation in the two directions is 

A = 2a/l = y/2a/b. (5) 

In our earlier paper Q, we have shown that the scat- 
tered flow from a sphere in a slit pore saturates ex- 
ponentially to an asymptotic Hele-Shaw velocity field 
at large distances from the sphere. In the asymptotic 
regime, the scattered flow thus assumes the form of a 
two-dimensional pressure-driven flow, with a harmonic 
pressure distribution which is independent of the trans- 
verse coordinate z. 

There are two ways in which the flow scattered from 
the small particles contributes to the hydrodynamic nor- 
mal force acting on the large particle. First, there is the 
near-field behavior of the scattered field, which produces 
stresses confined to a localized region of size defined by 
the length h ~ a. Second, there is also the far-field lu- 
brication pressure. 

To estimate the hydrodynamic force resulting from the 
near-field interactions, we note that there are N ~ Cl 2 a 
small particles in the lubrication region, and each particle 
yields an 0(r/aV p ) force contribution. We thus find that 
the cumulative effect of the near- field stresses on the total 
resistance force acting on the large particle is of the order 
0(4>7ibV p ), where 4> is the volume fraction of the small 
particles. This force is thus smaller by the order 0{a/b) 
than the lubrication resistance force 

F = 6nriV p b 2 /ho (6) 

resulting from the presence of the suspending fluid Q . 

The second small-particle contribution to the resis- 
tance force originates from the cumulative effect of the 
far-field pressure associated with the flow scattered from 
the particles. The far-field pressure decays slowly at large 
distances, and thus it acts over the whole lubrication 
area. As shown below, the far-field pressure produces an 
0(4>Fq) contribution to the resistance force. This contri- 
bution is thus the dominant correction to the hydrody- 
namic resistance. 

To find the far-field contribution of the scattered flow, 
we rescale the lateral coordinates s (parallel to the wall) 
and transverse coordinate z (across the gap) by the cor- 
responding lengthscales I and 2a, 

s — > s/l, z — » z/2a. (7) 

Accordingly, the Stokes equations ||2J| are transformed 
into the rescaled form, 

I dP X] fd 2 u z - 2 \ „ 
-2-a^ + i?{-0^ + XV H= F - (8a) 
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2a 



4a 2 



+ A 2 V V 



(8b) 



i=a 



2i , 



(11c) 



du z 
dz 



0, 



(8c) 



where subscripts s and z denote the directions parallel 
and normal to the wall, respectively. 

From now on, all coordinates and position vectors 
will be normalized by the corresponding characteristic 
lengths. To be consistent with this non-dimensional for- 
mulation, we introduce the following dimensionless pa- 
rameters: 



and 



&/2a=l/A 2 , h a = h /2a, 



5i = 5i/2a. 



h/2a, (9a) 



(9b) 



The subscript i in Eq. (|9bj) can be cither a (the small 
spheres) or b (the large sphere) or w (the wall) to in- 
dicate the excluded-volume range due to corresponding 
repulsive potentials. For simplicity, in our calculations, 
we take 



5 a = 5 W = 5b = 5/2 



(10) 



and use 5 as the only parameter which accounts for the 
excluded volume. The behavior of more general systems 
is similar to the one where Eq. (fTU)) is assumed. 

The rescaled Stokes equations (|5J| can be solved asymp- 
totically by scaling the velocity and the pressure fields 
with proper dimensional quantities and expanding these 
fields in the small parameter A, 



P = 



2a A 2 



V 1 



t£a^ 



(11a) 



(lib) 



i=0 



The different leading-order powers of A in the expansion 
of P, u s and u z stem from the different orders of various 
terms in Eqs. ©. By inserting expansion (fTTj) into ((SJ) 
and collecting terms with the same power of A, a hier- 
archy of equations corresponding to each order of A can 
be obtained. We are particularly interested in the zeroth 
order terms which satisfy the following leading-order lu- 
brication equations: 



^ = 0. 

dz 



- V S P 

du z0 
dz 



d 2 u M 
dz 2 
V s • u s0 



0. 



(12a) 



(12b) 
(12c) 



where 



. _ 4a 2 A 



(13) 



The leading-order hydrodynamic effect is obtained by 
solving the lubrication equations (fT2|) . 



B. Leading order pressure and friction 



The momentum equation in the z-direction (|12a[) im- 
plies that Po is independent of z. Using this fact along 
with the momentum equation in the lateral direction 
(|12bp and the no-slip boundary conditions at the solid 
surfaces, u s q can be determined in terms of the pressure 
field P , 



u s0 = \z(z - h)V s P + [ Z [ i s dz"dz' ~ Z T \\ i s dz"dz'. 
2 Jo Jo hj J 



(14) 



r 



Then, the velocity perpendicular to the wall can be eval- We combine Eqs. (fl"4"| - (fll)|) to obtain the lubrication 
uated by integrating the continuity equation (|12c|) equation for the zeroth order pressure field, 



u z0 



_ V.-u.,,/,'. (15) V s .(^V s P + b s ) = -12. 

The source term b s in Eq. (| 1 T[) is given by 
and using the boundary condition at the surface of the 

large particle b s (p) = B f s (r), 



Uzo\z=h = -!• 



(16) where 



(17) 
(18) 
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Bf, = -12 




i s dz"dz'~j-J J i s dz"dz'jdz. 



(19) 



In the absence of smaller species (i.e., when f s =0), 
well known lubrication equation for the pressure field in 
a slowly varying gap can be recovered from Eq. p7[) . 



V s • (h 3 V s P L ) 



42. 



(20) 



The presence of the smaller spheres produces a pressure 
correction 



Po 



Po-Pl 



(21) 



to the non-dimensional lubrication pressure Pl . The per- 
turbation pressure field po satisfies the lubrication equa- 
tion 



(22) 



Vs-(h 3 V s p + b s ) = 1 



which is obtained by combining Eqs. (fT7j) . (|20|) and (f2"Tj) . 

The function b s in the above equation plays a key role 
in our theory. It can be interpreted as the dipolar source 
density for the far-field lubrication pressure produced by 
the induced- force distribution F s . 

The correction in friction of the large particle due to 
the presence of the small spheres has contributions from 
both perturbation pressure po and viscous stresses cor- 
responding to the leading order velocity field Uo. How- 
ever, comparing the orders of A in Eq. (TTTJ we conclude 
that the second contribution is negligible compared to the 
first one. Thus, the leading order correction in A to the 
normal hydrodynamic force acting on the large particle 
in a solution of small particles stems from the leading- 
order pressure correction p . As a result, the correction 
in hydrodynamic force due to the presence of the smaller 
species can be expressed as 



2a 



77- / Pod 2 f 



(23) 



where po is assumed to vanish far away from the contact 
point. 



IV. DIPOLAR SOURCE DENSITY b s 



To evaluate the resistance force Fq from Eq. (f2"3"|) , we 
first need to calculate the lubrication pressure field po- 
The pressure can be determined from the lubrication 
equation (|22|) . provided that the dipolar source term (fl8|) 
is known. In this section, we focus on the key step where 
we compute the source term b s averaged over the posi- 
tions of small particles. 



A. Force singularities 

The induced force density at the surface of the z-th 
sphere is equivalent to a point-force singularity distribu- 
tion Fi at the sphere center. This description implies 
that the induced force density is considered at the sur- 
face of an infinitely small hypothetical sphere situated at 
the center of the i-th particle. However, the strength of 
the force density is modified in such a way that only the 
flow field inside the particle is altered, and the outside 
flow remains unchanged. Accordingly, the force density 
for an arbitrary particle i is 



4a 2 A 



^ ^ flmcr^lma ) j 



(24) 



where denotes the position of the center of the parti- 
cle. The coefficients fi ma represent the non-dimensional 
strength of the force singularity and are referred to as the 
multipolar moments. The force singularities Si ma can be 
expressed in terms of singular spherical basis solutions of 
the Stokes equations, 



V 2 v- 



= 5 



I ma • 



(25) 



The singular spherical basis solutions v JmCT and pressure 
solutions Pi ma are defined in [^, [1(3]. These fields sat- 
isfy the homogeneous Stokes equation at every point in 
the flow domain except at the particle center where they 
are singular and correspond to the point-force singularity 

By considering the cumulative effect from all the 
spheres, we can compute the statistically averaged nor- 
malized force distribution: 



F ( r ) = TTT — / ^( r o) y] fima{r )Si ma {r - r )d 3 r , 
4a z A na J ^— ' 



(26) 

where the volume fraction of small particles <fi is given by 

4-TT 

4> = Y aC - ( 27 ) 



The force distribution |26|) corresponds to the average 
induced-force density in the ensemble-averaged version 
of the Stokes equation ([2]) . Equation (|26[) is valid for an 
arbitrary distribution of small particles <f>. In our present 
application the volume fraction is uniform in the available 
space outside the excluded- volume shells. 
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B. Local parallel- wall geometry 

In order to evaluate the source term F(r), the multi- 
polar moments fi ma need to be computed. We obtain 
these multipolar force distributions by assuming that at 
the length scale a <C b, the region confined by the large 
particle and the wall can be locally approximated by the 
geometry of a channel bounded by two infinite parallel 
planar walls. 

For such a geometry, the multipolar moments fi ma can 
be evaluated with high accuracy either using a multiple- 
scattering technique [H[ or a Cartesian-representation 
algorithm developed in our recent papers 0, 

nana- 

The force multipoles fima(jo) are induced as a result of 
the interaction between the isolated freely moving small 
sphere situated at r and the incident horizontal lubrica- 
tion flow field uls created by motion of the large particle, 



d 2 u Ls 

dz 2 



= V S P L . 



(28) 



The lubrication pressure Pl satisfies Eq. ([20|) . so that 



From Eqs. ((28]) and J29J) , we find u Ls 



u Ls = 4A L t(l-t)e p) 



where 



3p 
Ah' 



(29) 



(30) 



(31) 



The amplitude Al depends on the radial position because 
both the volume flux and height of the gap are functions 
of p. 

At the length scale of the small particle diameter, the 
parabolic lubrication flow (|30|) can be treated as an ex- 
ternal flow coming from infinity in a parallel-wall chan- 
nel. It is convenient to choose a local coordinate sys- 
tem (a; ,y ,z) with the origin at the center of the small 
sphere and the x'-axis coinciding with the radial direc- 
tion e p . Hence, the direction of the incident parabolic 
flow is along the local x' coordinate. 

There are three immediate simplifications in the de- 
scription of fima for a freely suspended sphere in a 
parabolic flow. Firstly, the multipolar moments /i±io 
and /i±ii correspond to a Stokeslet and a rotlet, re- 
spectively, and therefore they vanish for a force-free and 
torque-free particle. Secondly, for a parabolic Poiseuille 
flow (like Ui„), the multipolar moments are non-zero only 
for m = ±1. Thirdly, our choice of the local coordinates 
implies that multipolar strength for m is real and is the 
same as the corresponding value for —m. In the follow- 
ing section we evaluate the relevant non-zero coefficients 
fi ma as a function of the gap width and position of the 
small particle, and we compute the pressure source term. 
C. The pressure source term 

In order to compute F§ from the integral ([23]). we need 
to determine po in Eq. (|22p . which involves the dipolar 
pressure source term h s . We find b s from f s by using Eq. 
(f!8|) . Combining Eqs. (f!3|) and (j26|) along with the proper 
non-dimensional scaling, we express f s in the following 
form: 



3& f 

f s (r) = — / 4>(r ) V" /; mCT (r )(I - e z e z ) • <5 im(T (r - r )d 3 r , 



(32) 



where I — e z e z is the projection tensor on x-y plane. 
Changing the order of the integrals in Eqs. (fT5|) and (p{2"|) 
we obtain b s : 

b s(p) = — [ 4>{r ) Y] iw(i*o)bi TO(T (p - p ; z )d 3 r , 
ira J f— 1 

Ima 

(33) 

where 

b;„ lCT (p - p Q ;z ) = (I - e z e z ) ■ B S hna (r - r ). (34) 

It can be shown that b; mcr is non-zero only when the 
following condition is satisfied Q 

I + cr - |m| < 2. (35) 

Due to the symmetry of the problem (a parabolic flow in 
a slit pore), only the terms with m = ±1 are relevant for 
our present analysis. 



In Appendix [A] we list the expressions for Si mrT with 
m = ±1 and derive the corresponding formulas for the 
source terms b; mcr , which are obtained by applying the 
operator B directly to 8i ma . Our analysis provides a 
simple derivation of the expressions for the far-field flow 
resulting from a force singularity in a slit pore. The pre- 
viously known formulas for the far-field behavior of a 
Stokeslet [l]| and rotlet [li| in a slit pore are special 
cases of our more general expressions. 

Considering Eq. (|3"5|) and m = ±1 for parabolic flows 
in a slit pore, we find that only bi±io, bi±n, bi±i2, 
b2±i0) t>2±ii, and b3±io are the nonzero contributions 
in Eq. (|33|) . Moreover, as pointed out in Sec. IIVB1 we 
have further simplifications: /i±io = /i±n = (freely 
suspended particle) and fi m a = fi -ma (choice of local 
coordinates). Hence, only four different coefficients fima- 
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(/ii2, /210, /211, /310) and eight different source terms 
b/ mCT (bi±i2, b 2 ±i , b 2 ±xii b 3± i ) are involved in the 
expression for the pressure source h s . In Appendix [Al we 
show that these source terms are given by the expressions 

b U a(r,r ) = b ; *_ lCT (r,r ) = d; CT (r )(e p + ie^,)8{p - p ), 

where the asterisk denotes the complex conjugate, 
6(p — po) is the Dirac delta function and 



dm 



d*2o = 6\/ — (h - 2z ), 



'hi 



^30 




-12 



47T 

4725' 



(37a) 



(37b) 



(37c) 



(37d) 



In the above relations, h denotes the local gap width, and 
e p and represent the unit vectors along the radial and 
azimuthal directions (note that these vectors correspond 
to the basis vectors e x > and e y i in the local coordinate 



system x' , y' , z' centered on a small particle). 

We combine the results (|3"rJj) and IpTF]) with Eq. 
find 



to 



b.(p) 



12/i 3 t 



(f>(p, zo)d(p, z )dz , 



(38) 



where d is the dipolar strength of the far-field flow pro- 
duced by a small particle at the position p, zq, 



d{p,z ) = 



2irah 3 



/,fn<r(P, Zo)dla{p, Z ). 



(39) 



This function is proportional to flow amplitude Al be- 
cause fn a varies linearly with A^. Moreover, d{p, zq) de- 
pends on the radial coordinate p only through the local 
gap width h(p) according to Eqs. (|37|) and the assump- 
tion that the induced-force multipolar moments fn a are 
evaluated in the local parallel-wall-channel approxima- 
tion. 
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FIG. 2: Normalized induced-dipole strength for two- 
dimensional Hele-Shaw pressure field as a function of position 
of the small-particle center from the wall for h = 2.0 (solid 
line), h = 5.0 (long dashed line), h = fO.O (dash-dot line) and 
h — 20.0 (short dash line). The dotted line represents the 
exact result for a very small particle in a parabolic flow when 
the wall is situated far away from the sphere. 



the multiple-scattering method lllll , but the Cartesian- 
representation technique d, [n], yields equivalent re- 
sults. In Fig. [2] we show only one half of the channel be- 
cause d is an even function of the particle position with 
respect to the channel mid-plane. 

We find that the maximum contribution in d comes 
from the multipole corresponding to the induced stresslet 
(I = 2, a = 1). When the particle is close to the wall, the 
induced stresslet increases to a large value because the 
slow inverse-logarithmic decay of the particle mobility 
creates a large hydrodynamic stress in the lubrication 
region. However, the dipolar strength d always remains 
finite, because for extremely small gaps the particle-wall 
relative velocity ultimately decreases to zero. 

There is no induced stresslet when the particle is at 
the center of the channel. For this position, only higher 
order multipoles contribute in d and therefore d has a 
minimum at the channel center. 

For large gap widths and intermediate positions of the 
particle where it is neither close to the wall nor near the 
center of the channel, the curves for different h match 
with the curve defined by 



D. Numerical results for the pressure source 

In Fig. [51 we present d/A^ as a function of the dis- 
tance of the center of the small particle from the lower 
wall for different gap widths where the position is normal- 
ized by the gap width. In our calculations the results for 
the force multipolar moments fima were obtained using 



d = b^-{l-2z/h) 2 . (40) 

This can be predicted by calculating the strength of the 
induced stresslet due to hydrodynamic interactions be- 
tween a parabolic flow and a sphere in an infinitely wide 
gap. 
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FIG. 3: Normalized integrated dipolar strength for two- 
dimensional Hele-Shaw pressure field as a function of normal- 
ized local gap for 8 = 0.1 (solid line), 8 = 0.2 (long dashed 
line), 8 = 0.5 (dash-dot line), 8 — 1.0 (short dash line) and 
8 = 2.0 (dotted line). 



FIG. 4: Correction to the lubrication pressure in the near- 
contact region due to the presence of small particles, nor- 
malized by Q-Krjbvp /a 2 , vs. the ratio between the local gap 
width and the small-particle radius, for different values of 8 
(indicated by the same line types as in Fig|3} ■ 



V. EVALUATION OF THE EFFECTIVE 
RESISTANCE FORCE F§ 

In this section, we focus on the evaluation of the cor- 
rection Fq to the hydrodynamic resistance coefficient for 
the large particle moving near the wall in a solution of 
small particles. To this end we first determine the pres- 
sure field for different ranges of the excluded shell 5. The 
resistance coefficient is then evaluated by integrating the 
pressure. 



We refer to D as the integrated dipolar strength. 

In Fig.[3j we show the magnitude of D as a function of 
the ratio between the gap width and the diameter of the 
small particle. Different curves correspond to different 
values of 5. When h/ (2a) < 1 + 25 the D is zero because 
there are no small particles in the gap. The magnitude of 
D however, increases rapidly to a maximum value when 



A. Evaluation of the pressure field 

For a given d, we can determine the resistance force F§ 
by combining Eqs. (|2"2"|) . (j2li|) and ([55]) . Here we perform 
explicit calculations for a model system of spheres inter- 
acting via the excluded-shell potential described in Sec. 
Ill Al Since in the short-time (high-frequency) limit the 
ensemble average in Eq. (I38p is taken over the equilib- 
rium particle distribution, the volume fraction of small 
particles is constant in the accessible region, 

(f>(r) = <f> when 1/2 + 5 < z < h - 1/2 - 6 (41) 

and otherwise (f> = 0. 

Using this model and Eq. we find 

b s (p) = 12 cj)A L h 3 De p , (42) 

where 

D= / l - S -'fSpV (43 ) 

J±+S A L 
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FIG. 5: Small-particle correction to the friction force acting 
on the large particle, normalized by Fq — 6nr]b 2 /ho, vs. the 
ratio between the minimum gap width at p — and the small 
particle radius, for different values of 5 (indicated by the same 
line types as in Fig0J. 
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the gap is wide enough to accommodate the particles. 
After reaching the maximum, the curves decay with the 
increasing gap because of the decreasing strength of the 
induced force multipoles. 

When a parabolic flow interacts with a sphere in a large 
gap, the multipolar moments can be calculated without 
considering the wall effects. In that case, the form (|4"0")) 
of the induced stresslet in the local shear flow yields 



D 



5 



(44) 



-2 in 



Figure [3] shows that the curves approach a slope of - 
the log — log plot for large h. 

Because of radial symmetry, all associated quantities 
in Eq. (|42|) like Al , h and D are functions of the radius p 
only. Hence, b s can be expressed in terms of the gradient 
of a scalar 



b s (p) = h 3 VS( P ), 



(45) 



where 



dp 



Therefore, by substituting Eq. 
the perturbative pressure po: 

Po = S. 



\2<j>A L D. (46) 
3 in Eq. QTZf), we find 



(47) 



As a result, the radial derivative of po can be obtained as 
a radial function by combining Eqs. (|4"6"| and (|47p along 
with the expression (|3 1 1) for amplitude Al of the incident 
parabolic lubrication flow 



dpo 
dp 



pD 
h 



(48) 



We determine po by integrating Eq. 

A further simplification can be achieved by noting that 
in the lubrication region the surface of the large particle 
is approximately given by 



Hp) = h o + Y- 



(49) 



Assuming (|4*9")l , the pressure can be expressed as an ex- 
plicit function of the local gap width: 



dpo 
dh 



-90- 



D 



(50) 



In Fig. [U the perturbation pressure field is plotted as 
a function of the gap width for different values of 8. The 
plateau region of the curves corresponds to D = in 
the region depleted of small spheres because of the ge- 
ometrical constraints. With the increasing range of the 
excluded- volume potential 6, the plateau region becomes 
wider, and a smaller number of particles is accommo- 
dated in the lubrication region. As a result p decreases 
with increasing 5. 



For a large gap, we find from Eq. (|44f that 



Po 



15(j> 
2h?' 



(51) 



which agrees with the result obtained using Einstein's 
formula for the effective viscosity. This asymptotic value 
is independent of the range of the excluded shell, pro- 
vided that 5 <C ft-o- This behavior is seen in Fig. 2] where 
for a large gap all the curves coincide with an asymptotic 
line of a slope of —2 in the log-log plot. 



B. Friction correction 

We now use Eq. (|2"5)) to obtain the correction to the 
normal friction coefficient of the large particle due to the 
presence of the smaller species. Owing to the axial sym- 
metry, the integral in Eq. (f2"5)) is reduced to the following 
form with the help of Eq. (|49|) : 



F c 



TTIjb 2 



(po-p^)dh, 



(52) 



ho 



The force F§, evaluated from Eq. (|52"|) . is presented 
in Fig. [5] as a function of the minimum separation hg 
between the wall and the surface of the large particle. 
The results are normalized by the product of the volume 
fraction 0, and the lubrication force ^ acting on the 
large particle in the absence of the small particles. The 
normalized correction to the friction force decreases with 
increasing 6 because fewer spheres are accommodated in 
the lubrication zone. For a larger gap, the ratio F§/Fq 
approaches a value of 2.50 for all values of 5. This is 
evident in Fig. [5] and consistent with Einstein's formula 
for the effective viscosity of the dilute colloidal solutions 
0. 



VI. CONCLUSIONS 

This article presents a new lubrication theory for ana- 
lyzing the near-wall dynamics of a large colloidal particle 
in the presence of smaller particles. Our theory is applied 
to determine the effect of the smaller species on the hy- 
drodynamic friction of the large particle when it moves in 
the direction normal to the wall. This effect is quantified 
in terms of the relative contribution F§ to short-time (or 
high-frequency) normal friction force acting on the large 
particle. 

We compute F§ as a function of the minimum separa- 
tion between the surface of a large particle and the wall. 
We model the screened electrostatic repulsion between 
the solid surfaces in terms of excluded volumes and vary 
the range of the excluded-volume potential to find its ef- 
fect on the normal friction. Our results show a decrease 
in friction with an increase in the excluded volume as the 
number of smaller particles decreases in the near-contact 
region. 
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The key simplification in our formulation stems from 
the fact that the predominant hydrodynamic effect of the 
small particles is due to the far-field form of the scattered 
flow from these particles. Our lubrication theory can be 
applied to evaluate this far-field velocity (which is ba- 
sically a pressure-driven Hele-Shaw flow in a slit pore). 
Our analysis provides a concise derivation of the general 
far-field solution for the flow produced by an arbitrary 
force singularity between two parallel walls. This deriva- 
tion is much simpler and more general than the earlier 
analyses of the Stokeslet [l3| and rotlet [l4[ flow in a slit 
pore geometry. 

We have verified our theory by comparing our calcu- 
lations with known limiting results. In particular, when 
the size of the smaller species is much smaller than the 
separation between the large particle and the wall, the 
calculated lubrication pressure field and normal friction 
correction agree with the expressions obtained by using 
Einstein's effective viscosity formulation Q. 

Our analysis can be extended to describe a more gen- 
eral near-contact dynamics of bidisperse suspensions. For 
example, our lubrication theory can also be applied to 
cases where non-spherical particles and non-planar walls 
are involved. In the future, we will focus on such a gen- 
eralization of the present analysis. We will also describe 
the effect of Brownian motion on the system dynamics 

El. 
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APPENDIX A: INTEGRAL DERIVATION OF 
FAR-FIELD FLOW IN A SLIT PORE 



2tt 8 2 8 



3 dz 2 ^ 



<4o(r) 



7T 8 2 8 



30 dz 2 6 ' 



4tt d 2 5 



4725 dz 2 



(Ale) 



(Aid) 



(Ale) 



(Alf) 



where 8(r) is the Dirac-delta function and the superscript 
|| denotes the projection of the vector on x-y plane. The 
vector e is defined as 



6 G_ T 



ie„ 



(A2) 



(except for normalization eis identical to e+i defined in 
[la]). For m = — 1, the source terms <5; mcr are the com- 
plex conjugates of 8n G . Therefore, we only focus on the 
force distribution described in Eq. (|A1|) . 

In the above expressions for 8i ma only 8, and the 
derivatives 88/ dz and d 2 S/dz 2 appear in addition to the 
combinatorial coefficients. Hence, in order to obtain b/ mcr 
by using Eqs. (|3"4")) and (|A1[) we derive the following re- 
lations by evaluating the integrals in the definition of 
operator B in Eq. (fT9|) : 



(A3a) 



B<5(r - r ) = 6z (h - z )S{p - p ) 

B-^-8(T-T ) = -6(h-2z )8(p-p ) (A3b) 
8z 

B^S(r-r ) = -128(p-p ). (A3c) 
8z z 



In this Appendix, we present the integral derivation of 
the far-field results described in section 3. As a corollary 
of our derivation, we recover the expressions for the well- 
known far-field flow solutions due to Stokeslet [l3[ and 
rotlet [li|. 

First, we list the horizontal component of 8i ma defined 
in Eq. (|24p . expressing the higher order multipoles as 
derivatives of Stokeslets. For our purpose, we only have 
to consider m = ±1 which corresponds to the parabolic 
flow. Moreover, the only multipolar forces 8i ma that are 
relevant are those which give non-zero b; m(T in Eq. (|34p . 
Hence, we take Eq. (|35[) into account and only analyze 
the cases where I + a < 3. 

For m = 1, the horizontal component of the relevant 
8ima is given below: 



*!ioW 



2tt 88 



3 dz 



(Ala) 



(Alb) 



By combining relations (|A1[) and (|A3|) with Eq. 
find 



we 



b lla (r,r ) = b;_ lCT (r,r ) = d l<J {r )e8{p - p ), (A4) 
where 

did = —2\f&nzQ(h — zq), (A5a) 



d n = -2Vfa(h~2z ), 



(A5b) 



and the remaining nonzero coefficients are given by Eq. 
(|37f . The relations (|A5a[) and (|A5b[) are associated with 
the far-field Hele-Shaw flow produced by a Stokeslet and 
a rotlet in parallel-plate geometry. Resulting expressions 
for 6no and 6m are in agreement with Liron-Mochon's 
[HI and Hackborn's [l4[ expressions for the far-field solu- 
tions whereas the other bi ma described in the article are 
equivalant to the coefficients derived in our earlier paper 
[7| by using a more complicated matrix representation. 
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